Supernarrow spectral peaks near a kinetic phase transition in a driven, nonlinear 

micromechanical oscillator 
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We measure the spectral densities of fluctuations of an underdamped nonlinear micromechanical 
oscillator. By applying a sufficiently large periodic excitation, two stable dynamical states are 
obtained within a particular range of driving frequency. White noise is injected into the excitation, 
allowing the system to overcome the activation barrier and switch between the two states. While the 
oscillator predominately resides in one of the two states for most excitation frequencies, a narrow 
range of frequencies exist where the occupations of the two states are approximately equal. At these 
frequencies, the oscillator undergoes a kinetic phase transition that resembles the phase transition 
of thermal equilibrium systems. We observe a supernarrow peak in the power spectral densities of 
fluctuations of the oscillator. This peak is centered at the excitation frequency and arises as a result 
of noise-induced transitions between the two dynamical states. 
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PACS numbers: 05.40.-a, 05.40.Ca, 05.45.-a, 89. 75. Da 



Periodically driven nonlinear systems often display 
multiple coexisting dynamical states under sufficiently 
strong driving fields. The presence of fluctuations en- 
ables these systems to occasionally overcome the activa- 
tion barrier in phase space, resulting in transitions be- 
tween the dynamical states 0. Noise- induced switching 
has been studied experimentally in a number of dynam- 
ical nonlinear systems that are far from equilibrium, in- 
cluding parametrically driven electrons in a Penning trap 
Q, radio frequency driven Josephson junctions and 
micro- and nano-mechanical oscillators 0, |H @ • When 
the noise is weak, the escape rate Ti out of state i (i = 1 
or 2) depends exponentially on the ratio of an activation 
barrier i?^ and the noise intensity I^: 



(1) 



Ri typically depends on system parameters such as the 
driving frequency and amplitude, as well as the shape 
of the power spectrum of the noise. The ratio of the 
populations of the two dynamical states is given by: 

-Ri)/In 



W1/W2 oc e^ 2 



(2) 



As a result of the exponential dependence of the popu- 
lation ratio on the difference in the activation barriers, 
the system will be found in either state 1 or state 2 with 
overwhelmingly large probability over most of the pa- 
rameter space y|, |jj . The occupations of the two states 
are comparable only over a very narrow range of parame- 
ters. This behavior bears close resemblance to systems in 
thermal equilibrium with two phases such as liquid and 
vapor. Such thermodynamic systems are usually in either 
one of the two phases and only at the phase transition 
will the two phases coexist. Even though driven, non- 
linear systems are in general far from equilibrium, theo- 
retical works predicted that a similar kinetic phase tran- 
sition would occur under the appropriate conditions Q. 
Like thermodynamic systems, fluctuations increase sig- 
nificantly when these non-equilibrium systems undergo 



phase transitions. A range of phenomena, including the 
appearance of a supernarrow peak in both the suscepti- 
bility and the spectral density of fluctuations [jj, |& li| , is 
expected to take place. 
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FIG. 1: (a) A cross-sectional schematic of the micromechan- 
ical torsional oscillator with electrical connections and mea- 
surement circuitry (not to scale), (b) Histogram of the oscilla- 
tion amplitude of the oscillator near a kinetic phase transition 
showing equal occupation of the two states, (c) Time trace 
of the normalized oscillation amplitude near a kinetic phase 
transition, demonstrating transitions between the two states. 

In this letter, we measure the spectral densities of 
fluctuations of an underdamped nonlinear micromechan- 
ical torsional oscillator near the kinetic phase transition 
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where the populations of the two attractors are compa- 
rable. The most prominent feature in the fluctuation 
spectrum is a narrow peak centered at the frequency of 
the periodic excitation. We demonstrate that this narrow 
peak is associated with noise- induced transitions between 
the two attractors. The width of the peak varies linearly 
with the transition rate and is more than a factor of 10 
smaller than the natural line width of the resonance in 
our experiment. Away from the kinetic phase transition, 
the intensity of the peak decreases exponentially. Apart 
from the narrow peak, we also observe smaller, much 
broader peaks in the spectrum that are associated with 
fluctuations within each attractor. These broad peaks 
are present for all driving frequencies within the hystere- 
sis loop and their dependence on the noise intensity is 
distinctly different from the narrow peaks at the kinetic 
phase transition. 

A micromechanical oscillator, consisting of a moveable 
polysilicon plate supported by two torsional rods, was 
used to take the measurements in our experiment. Two 
fixed electrodes are located beneath the movable polysil- 
icon plate, one on each side of the torsional rods. A 
voltage V applied to one of the electrodes generates an 
electrostatic torque that excites the torsional oscillations. 
The other electrode is used to capacitively detect the os- 
cillations. All measurements were taken at a temperature 
of 70 K and pressure < 10 -7 torr. The driving voltage 
V is a sum of a biasing dc voltage, a periodic ac voltage 
with frequency fd — ^ and random noise. The oscilla- 
tor, driven simultaneously by a periodic excitation and 
random noise, is well-described by the equation 



' + 2j9 + ulO + /36 3 = Esm(uj d t) + n{t), 



(3) 



where luq is the natural frequency for small oscillations, 
E is the effective amplitude of the periodic excitation, 
7 is the damping constant, n(t) is the effective noise in 
the excitation and /3 is the coefficient of the cubic non- 
linearity (an omitted term that is quadratic in 9 leads to 
renormalization of the nonlinear coefficients). Figure^, 
shows a schematic of this setup. A more thorough de- 
scription of the setup and derivation of the equation of 
motion can be found in Refs. H . 

In the absence of noise and at small driving torque, the 
response of the oscillator corresponds to that of a damped 
harmonic oscillator with a resonant frequency of 3286 Hz 
and Q w 8000. When the driving torque is increased 
beyond a critical value, the frequency response becomes 
hysteretic due to the cubic nonlinearity. Within a range 
of driving frequencies, two stable dynamical states coex- 
ist. In the absence of fluctuations, no transitions occur 
between the two states. When sufficient noise is injected 
into the driving voltage, the oscillator switches between 
the two attractors. Over most of the hysteresis loop, the 
activation barriers for escape from the two states 
are significantly different. The hollow and solid triangles 
in Fig. [2b shows the occupation of the high-amplitude 
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FIG. 2: (a) Normalized frequency response of the oscillator 
(circles) fitted to a damped oscillator with cubic nonlinearity 
(solid line), (b) Occupation of the two states versus frequency. 
At the kinetic phase transition (f p = 3284.8 Hz), the occu- 
pation of the high amplitude state (upright triangles) and the 
low amplitude state (inverted triangles) are comparable, (c) 
Dependence of the intensity of the supernarrow spectral peak 
on the driving frequency fd. The intensity falls off exponen- 
tially as the driving frequency is moved away from f p . 



and low-amplitude states respectively. On the low fre- 
quency side of the hysteresis loop, the occupation of the 
low amplitude state is considerably higher than the high 
amplitude state and the probability of finding the oscil- 
lator in the low amplitude state is close to unity. As the 
driving frequency decreases, the activation barrier R\ for 
switching out of the high amplitude states decreases. At 
the bifurcation frequency (/i = 3283.3 Hz), R\ goes to 
zero and the low-amplitude state becomes the only at- 
tractor. In an earlier experiment we showed that R\ 
depends on the detuning frequency (fd — fi) in the vicin- 
ity of the bifurcation point with a critical exponent of 
3/2 in agreement with theoretical predictions [lj. On the 
high frequency side of the hysteresis loop, a similar argu- 
ment applies except that the high-amplitude state is the 
stable attractor. 

While the oscillator is predominantly in one of the at- 
tractors over most of the hysteresis region, there exists 
a small range of frequencies where the occupation of the 
two attractors is of the same order of magnitude. Fig- 
ure ^ shows the oscillation amplitude as a function of 
time at a driving frequency of 3284.8 Hz and clearly il- 
lustrates the system switching between the two states. 
The relative occupation of the two states at this driving 
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frequency is deduced by calculating the area under the 
two peaks in the histogram of the oscillation amplitude 

(Fig. nt). 

Even though the nonlinear oscillator is a driven system 
that is far from equilibrium, the above behavior bears re- 
semblance to thermodynamical systems at phase transi- 
tions when two phases coexist. A kinetic phase transition 
was predicted to occur in non-equilibrium systems when 
the populations of the two dynamical states are equal 
0. An important feature associated with phase tran- 
sitions is the large fluctuations arising from transitions 
between the two states. We perform a systematic study 
of the fluctuations of the nonlinear oscillator at the ki- 
netic phase transition by examining the power spectral 
densities of the oscillator response. 

To resolve features in the spectral densities of fluctua- 
tions, it is necessary to have high resolution in frequency. 
The maximum number of consecutive data points that 
can be recorded with our instrument limits the frequency 
resolution in calculating the fast Fourier transform (FFT) 
of the angular displacement 6(t). Instead, we record the 
slowly varying envelope of the oscillations using a lock-in 
amplifier with a bandwidth of about 30 Hz. The re- 
sponse of the oscillator can be written as 

9(t) = A(t) cos(2tt f d t) + B(t) sm(2n f d t), (4) 

where A(t) and B{t) are the slowly varying amplitudes of 
oscillations in and out of phase with the driving torque 
at frequency fd- The spectral density of fluctuations is 
given by 

Q(/) = 4 EE [>*(*) -«*(*)) 

r t ( 5 ) 

x (A(t + t)+ iB(t + T))e- i2 ^ f ~ f ^ T ] , 

where TV is a normalization constant. 

In the absence of injected noise, oscillations occur only 
at the periodic driving frequency fd and the measured 
spectrum consists of a delta function centered at fd- 
When noise is added to the excitation, the spectral densi- 
ties of fluctuations become dramatically different. Figure 
13 shows the spectral density of fluctuations at 3 differ- 
ent periodic driving frequencies. To focus on fluctuations 
about the ensemble average response, the delta function 
peak at f d obtained with no injected noise is removed 
from the spectrum. In other words, the data point at f d 
is omitted for each panel in Fig. |31 Figure Oh shows the 
spectral density of fluctuations at the driving frequency 
fd ~ fp where the occupations of the two states are equal. 
The most prominent feature is a very sharp peak centered 
at the driving frequency. The width of this peak is a fac- 
tor of 10 smaller than the natural width of the resonance 
peak. This sharp peak is predicted to arise Q due to 
fluctuation-induced transition between the two dynami- 
cal states. Figures |2K and [3; shows the spectral density 



at two other driving frequencies that are comparatively 
far away from f p . As the periodic driving frequency is 
changed so that the oscillator moves away from the phase 
transition point, the sharp peak shifts accordingly to re- 
main centered at the driving frequency. The area under 
the peak, however, drops significantly. Figure |2t plots 
the area under the narrow peak as a function of periodic 
driving frequency, clearly demonstrating that the inten- 
sity of the supernarrow peak attains maximum at the 
kinetic phase transition and decreases exponentially as 
the occupation of one of the states exceeds the other and 
transitions between the states become less frequent. 
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FIG. 3: Power spectral density of fluctuations for three dif- 
ferent drive frequencies f d : (a) 3284 Hz, (b) 3284.7 Hz, (c) 
3285.1 Hz. Notice that the vertical scale of (b) is more than 
20 times larger than (a) and (c). In (d), (e) and (f) the axes 
have been rescaled to reveal the smaller and broader peaks in 
the fluctuation spectrums of (a), (b) and (c) respectively. 

Figure^ shows the behavior of the supernarrow peak 
with different noise intensities when the periodic driving 
frequency is held constant to maintain the oscillator at 
the kinetic phase transition. The noise intensity differs by 
a factor of 4 between solid and hollow circles. Both sets 
of data are fitted well by Lorentzians d, . As the noise 
increases, the peak width increases and the peak height 
decreases. We found that the area under the peaks re- 
mains about constant, changing by less than 10% when 
the noise intensity changes by more than a factor of four 
(left inset of Fig. The right inset of Fig. Et> plots 

the peak width as a function of the sum of the transi- 
tion rates Ti out of each state, where the transition rates 
are determined directly from the residence time between 
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transitions. The linear dependence of the peak width on 
the transition rate and the exponential decrease in the 
area of the peak away from the kinetic phase transition 
clearly identify the supernarrow peak with noise-induced 
transitions between the attractors. 
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FIG. 4: (a) The supernarrow peak for two intensities of in- 
jected noise at the kinetic phase transition. The injected noise 
intensity for the hollow circles is 4 times greater than the solid 
circles. The dotted and solid lines are Lorentzian fits to the 
data. Left inset: Dependence of the area of supernarrow peak 
on injected noise intensity. Right inset: The peak width vs. 
transition rate, (b) The broad, small peak at two injected 
noise intensities that differ by a factor of 2, at fd = 3284 Hz. 
Inset: Linear dependence of the peak area on injected noise 
intensity. 

In addition to the supernarrow peak, there are other 
peaks in the fluctuation spectrum that are weaker and 
much broader than the supernarrow peak at the kinetic 
phase transition. These peaks, unlike the supernarrow 
peak, are present for all excitation frequencies within 
the hysteresis loop, as shown in Figs. and |3f when 
the oscillator is far from the kinetic phase transition and 
in Fig. [3^ when the oscillator is at the kinetic phase 
transition. These smaller peaks represent characteristic 
frequencies of fluctuations about each dynamical state 
when the fluctuations are not strong enough to induce 
a transition over the activation barrier. In contrast to 
the supernarrow peak, the shape of these small peaks do 
not change considerably as the noise intensity increases 



(Fig. 0Jd) and their area varies proportionally with noise 
intensity (inset of Fig. QJd). 

Systems in thermal equilibrium, such as a Brownian 
particle fluctuating in symmetric double- wells potentials, 
also exhibit narrow peaks in their fluctuation spectrum 
|l l| . The broadening of these peaks with increasing noise 
intensity leads to the well known phenomenon of stochas- 
tic resonance [L^. We emphasize here the similarities 
and differences of our nonlinear oscillator to these sys- 
tems in thermal equilibrium. First, our oscillator is far 
from equilibrium [l3l llil Hjj and is not characterized by 
free energy. It is bistable only if it is driven strongly 
enough and is monostable otherwise. Second, at the ki- 
netic phase transition, the sharp spectral peak is centered 
at the driving frequency instead of zero frequency. The 
study of such critical kinetic phenomena in periodically 
driven nonlinear micro- and nano-mechanical oscillators 
could open new opportunities in tunable narrow band 
filtering and detection. 

We thank M. I. Dykman and D. Ryvkinc for useful 
discussions. 
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